* Plot IRF
local labely : variable label $depvar

local plotmin = 0
local plotmax = $hmax - 1

cap gen QuartersPlot = Quarters - 1 

if ${cfbands} == 1 {
	local lc1 = "${lc_1cfband}"
	local lp1 = "${lp_1cfband}"

	twoway ///
		(connected b QuartersPlot, lcolor(blue) lpattern(solid) lwidth(medthick) msymbol(oh) mcolor(blue) ) ///
		(line u${alpha1} QuartersPlot, lcolor("`lc1'") lpattern("`lp1'") lwidth(thin)) ///
		(line d${alpha1} QuartersPlot, lcolor("`lc1'") lpattern("`lp1'") lwidth(thin)) ///
		(line Zero QuartersPlot, lcolor(gs13)), legend(off) aspect(0.8, place(center)) ///
		title("{bf:Response of `labely' to a 1% of GDP tax cut}", color(black) size(medium) ) ///
		xlabel(`plotmin'(1)`plotmax' ) ytick( `yrange' , angle(0))   ///
		ytitle("Percent deviation", size(small)) xtitle("Quarters", size(small)) ///
		plotregion(margin(zero) color(white)) ///
		graphregion(color(white) margin(small)) 

}
else if ${cfbands} == 2 {
	local lc21 = "${lc_2cfband_alpha1}"
	local lp21 = "${lp_2cfband_alpha1}"
	local lc22 = "${lc_2cfband_alpha2}"
	local lp22 = "${lp_2cfband_alpha2}"

	twoway ///
		(connected b QuartersPlot, lcolor(blue) lpattern(solid) lwidth(medthick) msymbol(oh) mcolor(blue) ) ///
		(line u${alpha1} QuartersPlot, lcolor("`lc21'") lpattern("`lp21'") lwidth(thin)) ///
		(line d${alpha1} QuartersPlot, lcolor("`lc21'") lpattern("`lp21'") lwidth(thin)) ///
		(line u${alpha2} QuartersPlot, lcolor("`lc22'") lpattern("`lp22'") lwidth(thin)) ///
		(line d${alpha2} QuartersPlot, lcolor("`lc22'") lpattern("`lp22'") lwidth(thin)) ///
		(line Zero QuartersPlot, lcolor(gs13)), legend(off) aspect(0.8, place(center)) ///
		title("{bf:Response of `labely' to a 1% of GDP tax cut}", color(black) size(medium) ) ///
		xlabel(`plotmin'(1)`plotmax' ) ytick( `yrange' , angle(0))   ///
		ytitle("Percent deviation", size(small)) xtitle("Quarters", size(small)) ///
		plotregion(margin(zero) color(white)) ///
		graphregion(color(white) margin(small)) 
}
